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Abstract 

The bound state problem for a fermion-antifermion system is considered taking into ac- 
count a finite decay width of the constituents. We propose an exactly solvable relativistic 
zero order equation similar to that of Barbieri and Remiddi, but including a constant 
width. We focus especially on the tt system for which we reconsider our recent calculation 
of the bound state corrections to the toponium width, which was performed in the narrow 
width approximation and needed the use of second order Bethe-Salpeter perturbation 
theory. We show that one obtains the same result already in first order BS perturbation 
theory if one uses our present approach. Furthermore the large cancellations of gauge de- 
pendent terms is demonstrated to be a consequence of a Ward identity. This cancellation 
mechanism is shown to be valid for general fermion-antifermion systems. 
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1. The possible importance of off-shell corrections to the decay width in the calculation 
of the cross section for top-antitop production by e+e^ lately has been the subject of 
detailed analysis fT], 0]. In a recent publication we addressed this problem in the 
context of Bethe-Salpeter(BS) perturbation theory starting from the Barbieri-Remiddi 
equation 0. We were able to show by explicit calculation that large cancellations took 
place and that the net result can be interpreted as a pure time dilatation effect. This had 
been conjectured qualitatively before in view of the known results from muonium |^. 

In this note we will improve our understanding of this result and derive a general 
theorem, valid for similar cases. From a theoretical point of view our calculation for topo- 
nium has been unsatisfactory since it was obtained in the narrow width approximation. 
As we will explain below this approximation is not strictly applicable to that system. We 
now include the "bulk" of the decay width already in the zero order equation to avoid 
this problem. This offers also the more technical advantage that we then need only to 
consider first order BS- perturbation theory in contrast to the calculation in the narrow 
width approximation where second order perturbation theory was needed. Furthermore 
we are able to show that the cancellation of the gauge dependent terms is a result of a 
Ward identity. This fact also considerably simplifies the actual calculation. 

2. The Bethe-Salpeter (BS) approach for weakly bound systems is based upon the 
BS-equation 

G = D + DKG (1) 
for the ti Green function G{P, k, k'), where D is the product of the two (full) propagators 
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and K represents the 2pi BS-kernel. All four point functions in (|l|) depend on the total 
momentum P = (Pq, 0); the incoming (outgoing) lines carry momentum P/2±k {P/2±k'). 

Perturbation theory starts from an equation similar to eq. (|^) with Gq, Dq and Kq 
chosen in such a way that the exact solution is known. If Dq consists of the nonrelativis- 
tic fermion-propagators only and if Kq is the Coulomb kernel the zero order equation 
corresponding to (|l|) (after integrating out ko and k'^) simply reduces to the Schrodinger 
equation. Another very convenient zero order equation which already includes the rela- 
tivistic free fermion propagators and still remains solvable has been found by Barbieri 
and Remiddi (BR) some time ago However, as is the case for the top quark - if the 
decay width of the constituents becomes comparable to the binding energy, perturbation 
theory runs into troubles because of the occurrence of terms (r/(Q;^m))" in graphs like 
the one in fig. 1. To circumvent these difficulties one has to include at least a part of the 
exact self energy function in the zero order equation [Q. 

The key point of the BS approach as generalized to the case of a particle with finite 



width is the approximation to the free fermion propagator Ek = y k"^ + rn? 

S^{±^ + k,) = [(±^ + A:)7-m-S(±:^ + fc)]-i^ (2) 

^ ^ + ^ (3) 

±^ + kQ-Ek + i\ ±f + A;o + ^fc-2^' 

which is surely valid for a slowly moving particle. This choice leaves the relativistic pro- 



jectors 



2Ep 



unaffected and thus introduces a minimal number of terms involving the width F in the 



calculation. By integrating the zero order BS-equation for the Feynman amplitudes ( 

Do = 5+ ® 5_) 

tGo = -Do + DoKoGo (5) 

with respect to the zero components of momentum under the assumption that Kq is 
instantaneous and with a kernel of the special form {a = 4 /Sag) 



—Ana 

^2 



with 



/i = 2Ek/{Ek + m) 



V = 2/JPo + 2Ek + tT 



±7o) 



so as to annihilate the second term in we can show that 



(6) 



with 



Go = i(27r)^5(fc - fc')I^o + [A+A+A'+7o] ® [7oA"A-A'" 



X 



2uj 



X 



-nil 



{2Tif5{k-k') Gc{E,k,k') 



2uj 



muiy 



2J 



h.12 



/2 



Eh 



(7) 



is the solution of eq. (|^). Primes denote the dependence on the corresponding momenta 
and Gc{E,k, k') is the nonrelativistic Coulomb Green function [0] evaluated at 

Am 



This solution shows that the width of a particle can also consistently be included for a 
relativistic zero order equation similar to the nonrelativistic case where the replacement 
E ^ E + iV m the Green function has been proposed first |^. 

The corresponding solution for the BS-wave functions |^ {un = Ep — M^y2 — ie) 

Xn{p, e) = loXniP, -e)7o = f ^"^2^ ^i^^ A-w '^^^- 

{PO - ^n) nP) VM^"^ 

is identical to the BR wave-function for stable quarks and belongs to the spectrum of 




bound states P^'^ = M^°) — = y 1 — — iV. It appears as the real residue of the 
complex pole. The small imaginary quantity ie determines the integration around that 
pole. 

In eqs. (^) 5* is a constant 4x4 matrix which represents the spin state of the particle- 
antiparticle system: 

: singlet 

(10) 

<iml^~ '■ triplet. 

(j) is simply the normalized solution of the Schrodinger equation in momentum space, de- 
pending on the usual quantum numbers (n, /, m); a±i, ao in (|10D describe the triplet states. 
In the following it will sometimes be sufficient to use the nonrelativistic approximations 
of eqs. (D {uJn ^Ojn = {p'^/m + E„)/2 - ie) 

xip,er = ^^^ms = ionp,-erio (n) 

Po-^n 

Near a bound state, whose pole now lies in the lower half plane at Pq = the 
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Green function of the unperturbed (BR) equation behaves as 

G„iP,kX)^^^^^^ (12) 
A generic perturbation H changes Gq into 

oo 

G = Y.Go{HGor. (13) 

Since G, Gq, etc. denote Feynman amphtudes, H in (|T3|) differs from such an amphtude 
by a minus sign. It is straightforward to show that 

H = -K + Ko + iD-^ - IDq^ (14) 

The l.h.s of eq. ([131) looks hke (|12D with P^^ replaced by Pn and with perturbed wave 
functions. 

Formally expanding both sides of eq. (p!3| ) near a certain pole P^°^ 



Go = 9m{P, - Pi'Y^' (15) 

m=0 

oo 



m=0 



and using go = x ^ X allows the determination of the shift of the position of the pole 
APn = Pn — Pn^^ to arbitrary order 

AF 

AP„ = AM„ - = {{ho)) + {{hog,ho)) + {{ho)){{h)) + 0{h^). (16) 

A shorthand for BS expectation values has been used which incorporates four-integrations 
of the relative momenta { J dk = {2tt)~^ J d^k): 

((a)) = / dkdk'xij{k)aii'jj'{k,k')xi'j'{k'), (17) 



As already in eq. (0) etc., in composite expressions integrations and summations over 
internal variables are understood. i,j etc. represent spinor as well as color indices. A 
similar formula as eq. (0) holds for the perturbed wave functions It is important to 
note that eq. (^) is not a power series in a^. 

3. The correction to the decay width is obtained from the imaginary part of the 
expectation values, receiving contributions from the graphs in fig.2 and fig. 3. In ref. P| the 
problem of the gauge independent calculation of bound state corrections to the toponium 
decay width had been addressed. It was shown that it is possible to obtain a gauge 
independent result which can be interpreted as a correction due to time dilatation. But 
since the top quark decay width is of the same order of magnitude as the binding energy it 
is obviously preferable to use the equations valid for decaying particles as described above. 
In our previous approach it was essential that O(a^) corrections arose from all three terms 
denoted in ( [T^ ) because the BR equation was used to zero order. Here we will show that 
within our present formalism we need only to consider first order BS perturbation theory 
e.g. the first term in ([16|). 

Corrections to the fermion propagator have to be included via the perturbation kernel 
(see fig. 2) 



Po=Mi°'^-ir' 



We emphasize that the perturbation kernel has to be evaluated at Pq = Mj^^ — iT. With 
equ. (|l|) and ^lpo=Af<°'-ir = \=o 

/if^ = -z[(E(p-yn)+z-7o)®(;^'-m) + (;^-m)®(S(p'+yn)-z-7o)](27r)^5(fc-fc') (19) 



where n = (1,0) and p = ^^^^ + k, —p' = ^^f^ — k' denote the four momentum of the quark 
and the antiquark, respectively. The direct product refers to the t ® t spinor space. The 
factors (J) — m) and (J)' — m) compensate the superfluous propagator on the line without 
S. We have neglected terms that will contribute to 0(a^r^/m). 

Since we are calculating higher order effects we may also drop the dependence on the 
mass of the bottom quark although it may be included in principle. For the electroweak 
theory we use the i?g-gauge with the gauge fixing Lagrangian ( M denotes the mass of 
the W-Boson) 

'-gf = 

because it eliminates mixed W — (f) propagators. The gauge parameter ^ will not be fixed in 
the following. We obtain for the imaginary parts of E in (p!9|) (s = sm9w, V± = (1 ±75)72) 
for p2 > and p^^ > M^: 

-ImEvK = (20) 



C,f = -^\^^W:-^J<P^\' 



where 



^(0) = mY-^-'-^) (21) 
, jp' - M^f l-^ M^l + f 



Now we observe that the terms in (O) which include S have already been calculated 



in except for the imaginary parts occurring in the trace 

T := ^tr[7oA-5U+7oP±M^5A-(jf^-m)] 



^ l^{Vo + uj)[m{l - ^) + {Vo - uj) - I-]. (22) 

and in 

p — m ^ 2mpo — p — ma — imV. (23) 

We will be able to neglect these imaginary parts in the following since they either give 
rise to real corrections ( which we are not interested in) or to imaginary ones from ReS, 
both of which are of O^al^^^f,). Of course, these corrections would be required a in com- 
plete calculation of the toponium decay width to this order in order to obtain a gauge 
independent result. This goes beyond the scope of the present paper. Furthermore, these 
corrections are state independent, expected to be very small and therefore of minor prac- 
tical interest, anyhow. As can be seen from ( ^2]) the axial part of S does not contribute 
within the expectation value. 

Therefore, the contribution from first order BS perturbation theory to the decay width 
to O(a^r) reads 

((/.«)) = ri + r3 (24) 

with 

Ti = -4lm{{-tJ:{p) ® {f - m){27v)*6{k - k'))) = (25) 

= To - 2T,al + ^[(2M^ + m^)^^^^^ + 2p(m^0](^^), 
Ts = -4lm{{^^{^'-m){27rY6{k-k'))) = -To{l + ^) (26) 

A factor 2 arises here and in the following, counting both self energy contributions of the 
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t and t line: 

In the sum Fi + the leading Fq terms cancel and the resulting correction is already of 
0(a|Fo). 

Now we turn to the calculation of the vertex correction (Fig. 3). Performing the color 
trace this gives rise to the kernel 

hf = tlmAo ® (28) 

In Q we calculated ImAo := limq^o,p2^m2 A(p, q) directly. Here we will show that it is 
possible to use a Ward identity to determine this correction. We observe that 

^o(P,9)|^^o = -gQCDT'^-^Eip) (29) 

holds at least at the one-loop level. In the graphs shown in fig. 3 to the required order it 
proves sufficient to evaluate the leading contribution by putting the fermions on the mass 
shell and applying the nonrelativistic wave functions ([TT]). Using the exphcit expression 
(pop for ImE we obtain 



loF (30) 



327rs2 

(1 + -—)(! + 3—) + 4p(^', 0- 



This result agrees with that of ref. and one obtains a correction to the decay width: 

Ana, 

(^)^- (31) 
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The fact that the gauge dependent terms in the sum Fi + cancel is thus traced back 
to the identity (pQl). The gauge independent contribution from first order BS-perturbation 
theory is now different from that of |Q (actually it gives now the net result) in accordance 
with the different zero order equation, used in our present context. 

To complete our calculation we finally check the contributions of second order pertur- 
bation theory. It can be shown that due to the fact that the first weak correction now is 
already of order 0{alr), this corrections do not contribute to the required order. Thus 
our present approach simplifies the calculation of bound state corrections to the decay 
width considerably. 

Summing up all contributions we get the result 

^^boundstate = Tl + + Ts = ""Tq-^ (32) 

in agreement with our previous calculation 0. 

4. The above considerations suggest that the result (P^ ) is only a special case of a more 
fundamental statement, with broader range of applicability, which we will now derive. 

Consider a fermion whose decay can be described by the imaginary part of a self energy 
function which will have the general form (in any covariant gauge in the relevant sector 
of the theory) 

E(p) = S5(/) +^Ev{p')+-f5^p{p')+^,^A{p'). (33) 

Within the expectation value (|T7|) the pseudoscalar and axial vector parts vanish. Fur- 
thermore we will drop all factors Q{p^ — fi"^), which is a valid approximation in all cases 
where a^m <^ m — yU. If this is not the case, either the whole calculation makes no sense 
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because then one would have to include the entire rung of Coulomb interactions of the 
remaining particle with the heavy decay products or the decay width becomes very small, 
or both. Let us therefore consider the correction to the decay width which results from 
the remaining parts of the self energy in an on shell renormalization scheme ( m is the 
pole mass of the fermion) . We expand S around that mass shell: 

= z(ImS5 +|^ImSy) + {J^'s + ^^y) {p^ - m^). (34) 



Terms of 0(a;^g^;.) are understood to be neglected (cf. remarks after eq. p3| ) ) and we 
denote 

Ex := Sx(m2), X = S,V 
d 



2^ 



From eqs. ( pQ]) and ( |3lD we can easily calculate the vertex correction, say, for the upper 
particle line: 

AoA+ = ^T"7o[zImSy + 2m{J:'s + mSy)]A+ (35) 
After some algebra we can write to the required order 

hf^ = (ImScj + mImSy + -70) - m) - iD^l^Q [zImSy + 2m{J:'s + mJ^'y)] 

(36) 

ho"^ = -^p[zImSy + 2m{T!g + mSy)]7o ® 70 



Since 

r = — 2MlmEM 



^ = — ImSo(m^) — mlmSy(m^) (37) 

p=(m,0) * V ; v\ / \ J 
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and (ziI'or=o ~ -^o,r=o)x = we obtain for the bound state correction to the decay with 
to Oia^a^eak) 

Ar^ow^tate = -4^r((A- 0if- m))) = -ro^ (38) 

Now it is clear from the preceding argument that it apphes not only to the toponium 
system but also to all other systems which can be described as weakly bound fermion- 
antifermion system with unstable components. Another example which is known for a 
long time is provided by the muonium system where the above result has also been 
obtained first by explicit calculation [Q. We conclude that the formalism developed here 
and especially the use of the identity ( ^91) simplifies the problem of the bound state 
correction to the decay width in a profound way. It is now possible to clearly isolate the 
underlying cancellation mechanism which automatically gives a gauge independent result 
which can be interpreted as time dilatation alone. However, because we may now also 
include the scalar part of the self energy function and since Kq can be any kernel, the 
result ( pSf ) is found to comprise and as special cases of a more general theorem: 
The leading bound state corrections for weakly bound systems of unstable fermions (with 
decays like t — >• 6 + W^,n~ — >• e~ + z/g + t'^) are always of the form (P5|). Among the 



consequences we especially not that (|58|) may be even safely applied to toponium even if 



one uses e.g. a renormalization group improved potential. 
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t -i(p'-m) p'= - P/2 + k 

Fig. 1: Multi-decay contribution 




- p/2 + k 

= - p/2 + k 



Fig. 2: Weak oflf-shell corrections ( D"* ) 




Fig. 3: Coulomb enhancement { -Iho*^ ) 
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